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We describe a unified approach for the determination of ac-Stark correction and Kramers- 
Heisenberg dispersion formula. In both cases the contribution from infinite intermediate summa- 
tion appearing in the expression for the corresponding matrix elements are evaluated exactly in the 
dipole approximation for the ground state of hydrogen atom using a variation of the Dalgarno-Lewis 
method. The analytical expressions obtained can be efficiently used for the numerical evaluation 
of matrix element for all values of incident photon energy and comparison is made with results 
obtained by different methods. 

PACS numbers: 32.80.Rm, 42. 50. Hz 



I. INTRODUCTION 

In the presence of intense laser field, atomic or molec- 
ular levels are displaced (shifted or broadened) and these 
stimulated radiative corrections are alternatively known 
as ac-Stark shift, dynamic Stark effect or light shift 0. 
Being a fundamental process ac-stark effect is very well 
studied theoretically as well as experimentally in wide 
area of atomic and molecular physics. This effect is also 
very important in the area of laser trapping and cooling 
and is a basic mechanism behind many nonlinear opti- 
cal effects |2J . It is a promising tool for the field of op- 
tical communications (Optical switching methods) and 
plays a very decisive role in many current research ar- 
eas S 0, II SSSIIillllEEH In this work 
we have described an efficient method to obtain a closed 
form analytical expression for two of the basic process 
viz. ac stark shift and elastic scattering cross section 
of photons for the atomic hydrogen in the ground state. 
Our aim is to get an analytical expression for dipole dy- 
namical polarizability (uj) from which one can cal- 
culate the level shift SE^ and level width T^ 2 \ The 
hydrogen atom being the simplest of all atomic systems, 
plays a very special role in this respect because it has 
closed-form wavefunctions. Since these wave functions 
are analytically known, both the dipole dynamic polar- 
izability and Kramer's-Heisenberg dispersion formula for 
hydrogen atom can be written in a closed form. 

The level shift for atomic hydrogen was previously cal- 
culated by both per turbative and nonperturbative meth- 
ods [H [H [H E HI IH HI . In the per- 
turbative method, evaluation of the higher order matrix 
element is the key problem and it is achieved through 
Greens function formalism [l6l Il7t ITs| or b y so lving a 
system of second-order differential equations 20] . 

The Kramers-Heisenberg formula alternatively known 



as dispersion formula plays a very important role in the 
determination of scattering cross-sections |25j. Consid- 
ering its importance as an age-old relation there were 
only very few attempts to evaluate it analytically |26ll27| . 
There were some generalization of it and alternate forms 
are also available in the literature [2^ . I30L l3ll l32| . Re- 
cently a low-energy expansion of the Kramers-Heisenberg 
formula for Rayleigh scattering process in atomic hydro- 
gen was reported |33(. Similar to the ac Stark shift this 
also involves the evaluation of the second order matrix 
element. 

The exact calculation of higher order process in the 
perturbative formulation is nontrivial because of the 
presence of infinite summation over intermediate states 
in the expression for the higher order matrix element. 
In our approach the matrix element containing infinite 
summation over the whole hydrogenic spectrum (discrete 
and continuum) is performed by using an implicit sum- 
mation technique due to Dalgarno Lewis |34| . which re- 
duces the evaluation of the infinite summation to finding 
a solution of some inhomogeneous differential equation. 
The closed form expression which we have obtained is 
very simple and also very convenient for analytical con- 
tinuation. Thus with this method, wc can very easily 
obtain the relevant matrix element for radiation with en- 
ergy larger than ionization energy (above threshold ion- 
ization), while other methods need some kind of approx- 
imation like Pade approximation in the case of Coulomb 
Greens function formalism. 



II. RADIATIVE CORRECTION IN THE 
GROUND STATE 

The level shift depends on the intensity / and the fre- 
quency uj of the radiation and the complex second order 
shift AW is given by [II [II [11111 
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A^(cu) = 6E^ - iT^ = -± t W(uj) (1) 
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where the real part SE^ is the energy shift and the imag- 
inary part I^ 2 ) gives the level width and Iq = 7.016 x 10 16 
W/cm 2 is the characteristic atomic field strength inten- 
sity. Here (uj) is the dipole dynamic polarizability 
and in the atomic unit it can be written as 



*(p, q, A, r) = / ds e- r ^K(p, q, A, s) (11) 
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{E g -E n -u) 
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where E g is the atomic ground state energy and e is the 
polarization of the radiation. First term in the bracket 
is the absorption-emission term and the second is the 
emission-absorption term. The infinite summation over 
the complete set of intermediate state in equation @ can 
be performed exactly by defining a set of operators F and 
F such that 



e • r \g) = {FH a - H F + loF) \g) 



e*-r\g) = (FHo-HoF-ujF) \g) 



(3) 
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K(p, q, A, s) = 



K{p, q, A, s) 
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(12) 
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where A = yT — 2w and A = yT + 2cj. In this work we 
obtain the solutions with integer values of p and g. But 
in general for the purpose of analytic continuation it can 
be complex and A also can become complex depending on 
the value of the frequency u>. Now using equations i|l(J|) 
and l|ll|) in equation (JSJ the final form of the dipolar 
polarizability becomes 



where Hq = — V 2 /2 — 1/r is the unperturbed Hamil- 

7r is the 



tonian for atomic hydrogen and \g) = e 
ground state wave function in atomic units. With these 
definitions and the closure relation | n)(n \ = I the 
expression in equation for dipole dynamical polariz- 
ability will be reduced to 



T<»(u) = {g\e* ■ rF \ g) + {g \ e ■ rF \ g). 



(5) 



Thus the infinite summation over the intermediate states 
are reduced to the determination of the operator F and 
F. Using co-ordinate space representation and writing 
F = e ■ r/(r) and F = e* ■ r/(r), the equations |(3J| and 
Q become 

r ~T^f( r ) + (4 - 2r) — f(r) + {2ur - 2)/(r) - 2r (6) 

r $2^ r ) + ( 4 - 2 r)T-/M ~ (2wr + 2)/(r) = 2r (7) 
dr dr 

By the method of Laplace transform 36] for the solution 
of differential equation we can obtain the solutions to the 
above differential equations as 



and 



with 



/(r) = --^ 5 #(l l l,A,r) 
lo 2w J 



/(r) = 1 + ^5(1, l,A,r) 



(8) 



(9) 



*(p, g, A, r) = / ds e- r ^K{p, q, A, s) (10) 
x 



T ( 2 )( w ) = — / dre- 2r r 4 {$(l,l,A,r)-$(l,l,A,r)} 
3W 3 7 
o 

(14) 

The limit ui — > of the expression for r"'(a;) can be 
shown to approach the value corresponding to the static 
dipolar polarizability which is 9/2. 



III. ELASTIC SCATTERING OF PHOTONS 

The differential scattering cross section for low energy 
elastic scattering of photons with fre q uen cy u >, by atoms 
is given by the dispersion formula [25l l26l |27j 



a 2 ( e . e ') 2 |i-^H-^(-^)| 2 



where ao is the Bohr radius, M is the Kramers- 
Heisenberg matrix element, e and e' respectively are the 
initial and final polarization of photons and P(u>) in 
atomic units is given as |27j 
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Here p is the momentum operator and the summation 
is over the complete set of states including continuum 
states. It is useful to note from equation (|15|l that the 
differential cross section for coherent scattering of pho- 
tons is just the Thompson cross section modified by the 
dynamic polarizability. In a similar fashion we can con- 
sider Raman scattering where initial and final states are 
different. The analytical expressions for P(u>) were de- 
rived earlier using Schwartz and Teinman method j2(J 
and Coulomb Greens function (CGF) formalism p7j . 
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Using a slight variation of the formalism in the previ- 
ous section we can easily calculate P(ui). In this case the 
infinite summation over the intermediate state in equa- 
tion (|16f) is performed by defining a set of operators U 
and U such that, 



integrals it is very convenient to define the following 



p\g) = (VH -H V + uV)\g) 
p \g) = (lJH - #oU - ccU) \g) 

and the expression for M will become 
2 



(17) 
(18) 



M(cu) = 1 + - (fo | p • U | g) + {g | p • U | g)j (19) 

Now the evaluation of the infinite summation is reduced 
to the evaluation of the operators U and U. To obtain 
a similar expression as in the previous section, instead of 
going to momentum space representation [2a,|22l, we use 
the coordinate space representation of U and U. This 
is done by taking U = ru(r) and U = ru(r) and the 
equations (|17|l and (|18H become 

d 2 d 
r—ru(r) + (4 - 2r) — u{r) + (2ur - 2)u(r) = 2i (20) 



r—ru(r) + (4 - 2r) — u(r) - {2wr + 2)u{r) = 2i (21) 
dr dr 

These differential equations has the same form which ap- 
peared in the previous section and the solutions can be 
written as 



u(r) 



2uo 



2 $(l,l,A,r) 



and 



u{r) = ^$(l,l,A,r) 



(22) 



(23) 



Using this in the definition of U and U and substituting 
it in equation (|19|) we get a closed form expression for 
the Kramers- Heisenberg matrix element. It now takes 
the form 



I(p,q,X,n) = j dre ZT r n $(p,g, A,r) 
n 



(26) 



= r,\ 



1 - A 



i i 



( s + xf +i (s-xy 



I + XJ J " (l + s)™+ x 

A 



and 



I(p,q,X,n) = dre ZT r n $(p, q, X,r) 



(27) 



A + l 



(l + s )n+l 



With these, the Kramers-Heisenberg matrix element and 
the dipole dynamic polarizability respectively will take 
the following forms 

M(w) = 1-^2 A > 3 ) - 'C 1 ' l >*> 3 )} ( 28 ) 



r {2) (w) - ^ { J(l, 1, A, 4) - 7(1, 1, A, 4)} (29) 



This is the final form of the analytic expression for the 
matrix elements. We did not attempt to rewrite these 
in terms of well known special functions because the re- 
maining integrals in equations (|26|l and (|27|l can be done 
very efficiently by a simple numerical integration. 

The integrals in equation (|10fl is defined only for Re(g— 
1/A) > — 1. But in order to calculate the matrix elements 
in equations (|28|l and (|29() for all values of w, we require 
an analytical continuation of the definition of $(.) given 
in equation 11U|) . It can be achieved by performing a 
partial integration in the definition given in equation I|1U|) 
and this gives 

1 



MM = l-i /dre- 2 »-r 3 {$(l,l,A J r)-*(l,l,A,r)} $ ^' A < r ) = q + 1 _ i ^ 9 + 1, A) + r$(p, g + 1,A(3Q) 



(24) 

If the incident photons are unpolarizcd and the polariza- 
tion of the scattered photons are not observed the differ- 
ential scattering cross section will take the standard form 



dcr = r o\( l + cos 2 9)\M(uj)\ 2 dn 



IV. DISCUSSION AND CONCLUSION 



(25) 



The radial integrals in equations l|14|) and (|24|l can be 
done exactly |37| . For the numerical evaluation of various 



A 

!> + - )<!>(<<- + :i. A./-) 



This is a very important recurrence relation because l|3U|) 
can be used to analytically continue the definition of 
these function for Re(q — 1/A) > —2. We can also re- 
peat this procedure to analytically continue the definition 
to other values of Rc(q — 1/A). And for the numerical 
evaluation of various integrals as a function of radiation 
frequency w we have extensively used this relation. It 
is also very important to note the appearance of simple 
poles in the definition of $(.) when q = — 1 + 1/A. Since 
A = \/T— 2oj, this is due to the presence of intermediate 
resonances (when 2u> = 1 — 1/n 2 , where n = 2, 3, . . .) in 
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the expression for second order matrix element in the per- 
turbation theory. This also leads to a strong frequency 
dependence of ac stark effect and scattering cross section. 

Similarly for above one photon ionization threshold i.e. 
lu > 1/2, a simple analytic continuation makes A to be a 
purely imaginary number. Thus by this simple method, 
with the same analytic expression for the matrix element, 
we can numerically evaluate it for the whole physical 
range of photon frequency ui. Displayed in table [I] are 
t' 2 ) (uj) as a function of photon frequency in atomic units 
and are compared with values obtained by other meth- 
ods. It is purely real for lu < 1/2. For u> > 1/2 one 
photon ionization is possible and this makes the inter- 
mediate virtual state to lie in the continuum and t^ 2 -* (lu) 
becomes complex. 

The dipole dynamic polarizability for photon energy 
below and above one photon ionization threshold is given 
in tabled Our values are compared with values obtained 
by Arnous et al. |17| using Coulomb Greens function 
method, but with an additional factor of lu~ 2 in their 
result. It is very easy to see that this factor is missing 
from their result by taking the lu — > limit. In this 
limit our result approaches the dc polarizability, which is 



9/2. The values given in [l7| also approaches to the same 
limit provided their results are multiplied by an overall 
factor of lu~ 2 . Using the analytic continuation described 
in equation (|3(JI) we can calculate (lu) for lu very close 
to intermediate resonance also. Our result also agrees 
with the second order level shift and width reported by 
Pan et al. |lfi| . From the values given in tablesQand|nlit 
is useful to note the change in sign of the level shift when 
the photon frequency cross these resonance values. The 
values of Kramers-Heisenberg matrix element are given 
in tables [H] and II I II They are in good agreement with 
the results of Gavrila |27) . 
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TABLE I: Comparison of the values of dipole dynamical polarizability in atomic units for both below and above one photon 
ionization threshold lo = 0.5. In the limit to — > 0, t^- 2 \uj) approaches the dc polarizability, which is -4.5. 



oj(a.u) 


This work 


ABM " 


cj(a.u) 


This work 


ABM a 


This work 


ABM" 


0.001 


4.50003 


- 


0.6 


3.297 


- 


2.505 


- 


0.002 


4.50011 


- 


0.7 


2.493 


- 


1.408 


- 


0.02 


4.51066 


-4.51 


0.8 


1.915 


1.915 


0.850 


0.8506 


0.04 


4.5429 


-4.5431 


1.0 


1.205 


1.205 


0.362 


0.3627 


0.08 


4.6775 


-4.6776 


2.0 


0.275 


0.275 


0.023 


0.0239 


0.10 


4.7843 


-4.7843 


3.0 


0.117 




0.004 




0.20 


5.9416 


-5.9416 


4.0 


0.064 




0.001 




0.43 


0.2971 




5.0 


0.041 


0.041 


0.0005 


0.00057 


0.46 


3.9273 




6.0 


0.028 




0.00027 




0.465 


3.0867 




9.0 


0.012 




0.000049 




0.477 


1.2644 




10 


0.010081 


0.01008 


0.0000319 


0.00003 


0.478 


1.9330 














0.489 


0.6465 















"Values obtained using Coulomb Green function Il7l by Arnous 
ct al. with an additional u>~ 2 factor as explained in section ft VI 



TABLE II: Kramers-Heisenberg matrix element M in atomic units for lo below one photon ionization threshold. Comparison 
is made with values in Ref. 12711 . 





This 






This 






This 


oj(a.u) 


Work 


Ref. [27] 


w(a.u) 


Work 


Ref. [27] 


oj(a.u) 


Work 


0.002 


0.000018 




0.376 


77.8416 








0.02 


0.0018 


0.0018 


0.38 


15.3829 


15.3828 


0.481 


4.0681 


0.04 


0.0072 


0.0072 


0.4 


2.6916 


2.6916 


0.484 


0.6692 


0.06 


0.0165 


0.0165 


0.429 


0.0611 




0.485 


0.4490 


0.08 


0.0299 


0.0299 


0.43 


0.0549 


0.0549 


0.486 


16.087 


0.10 


0.0478 


0.0478 


0.44 


3.1503 


3.1503 






0.12 


0.0708 


0.0708 


0.444 


38.8927 




0.488 


1.2367 


0.14 


0.0999 


0.0999 








0.489 


0.1546 


0.16 


0.1361 


0.1361 


0.445 


32.2604 


32.2603 


0.49 


6.6498 


0.18 


0.1812 


0.1812 


0.453 


2.3124 


2.3124 


0.491 


1.2466 


0.20 


0.2376 


0.2376 


0.464 


0.2004 




0.492 


3.0681 


0.22 


0.3091 


0.3091 


0.465 


0.6674 


0.6674 


0.493 


1.2572 


0.24 


0.4016 


0.4016 


0.468 


8.1693 


8.1693 


0.494 


3.9947 


0.26 


0.5246 


0.5246 








0.496 


3.0200 


0.30 


0.9507 


0.9507 


0.469 


27.9814 


27.9814 


0.497 


3.1600 


0.32 


1.3752 


1.3752 


0.473 


1.97857 


1.9785 


0.497 


3.1600 


0.36 


5.3036 


5.3036 


0.477 


0.2876 




0.498 


0.7238 


0.37 


15.763 


15.763 


0.478 


0.4416 


0.4416 
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TABLE III: Real and Imaginary part of Kramers-Heisenberg matrix element for photons of energies above one photon ionization 
threshold (lo > 0.5) and comparison is made with values in Ref. [27) . 



Re M Im M 



oj(a.u) 


This work 


Ref. [271 


This work 


Ref. [27] 


0.6 


1.1872 


1.1872 


0.9018 


0.9018 


0.7 


1.22161 


1.2216 


0.6900 


0.6900 


0.8 


1.22612 


1.2261 


0.5444 


0.5444 


0.9 


1.21842 


1.2184 


0.4400 


0.4400 


1.0 


1.20598 


1.2059 


0.3627 


0.3627 


2.0 


1.10007 


1.10007 


0.0958 


0.0958 


3.0 


1.05696 


1.0569 


0.0421 


0.0421 


4.0 


1.03685 


1.0368 


0.0231 


0.0231 


5.0 


1.02589 


1.0258 


0.0144 


0.0144 


6.0 


1.01924 


1.0192 


0.00977 


0.00977 


7.0 


1.01489 


1.0148 


0.00699 


0.00699 


8.0 


1.01188 


1.0118 


0.00522 


0.00522 


9.0 


1.00971 


1.0097 


0.00403 


0.00403 


10 


1.0081 


1.0081 


0.00319 


0.00319 


20 


1.00236 


1.0023 


0.00066 


0.00066 


30 


1.00112 




0.000262 




40 


1.00066 




0.000133 




50 


1.00042 




0.000075 




90 


1.00014 




0.0000196 





